Control of the Bose-Einstein condensate by dissipation. Nonlinear Zeno effect 
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We show that controUed dissipation can be used as a tool for exploring fundamental phenomena 
and managing mesoscopic systems of cold atoms and Bose-Einstein condensates. Even the simplest 
boson- Josephson junction, that is, a Bose-Einstein condensate in a double-well trap, subjected to 
removal of atoms from one of the two potential minima allows one to observe such phenomena as the 
suppression of losses and the nonlinear Zeno effect. In such a system the controlled dissipation can 
be used to create desired macroscopic states and implement controlled switching among different 
quantum regimes. 
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I. INTRODUCTION 

Universal properties of condensed atomic gases as sys- 
tems hosting spectacular nonlinear phenomena like insta- 
bilities and collapse, solitons and shock waves, localized 
modes, vortices and self-trapping and delocalizing tran- 
sitions, nowadays are well known These phenomena 
constitute the properties of Hamiltonian systems, where 
dissipation is considered as an undesirable destructing 
factor. This role of dissipation can be inverted, if a 
system possesses an intrinsic mechanism balancing the 
losses and giving origin to stable dissipative structures 
. Alternatively, the dissipation can lead to a construc- 
tive effect when its action is limited in time, allowing one 
to generate diverse nonlinear excitations [3|, or when it 
has a nonlinear origin Q supporting localized patterns. 
These effects appear on the macroscopic scale. Study 
of the dissipative decay can reveal also the microscopic 
quantum properties of the atomic gases Q and even in- 
hibit the losses of atoms, by inducing strongcorrelations 
due to a large imaginary scattering length [g . 

The purpose of this article is to show that there are 
various dissipative regimes in systems of cold atoms and 
Bose-Einstein condensates (BECs) loaded in a multiwell 
trap and that removal of atoms can serve as a tool for 
exploring the fundamental quantum phenomena. In par- 
ticular, we find that the dissipation can inhibit losses and 
allows one to manage the so-called Macroscopic Quan- 
tum Self- Trapping state (below simply the self-trapping 
state) 0, [1| . Moreover we observe the inhibition of losses 
due to atomic interactions which can be termed as the 
macroscopic nonlinear Zeno effect. 

Specifically, we consider the case of a BEC in a double- 
well potential whose Hamiltonian dynamics is well un- 
derstood theoretically and has been a subject of fun- 
damental experiments and numerous potential applica- 
tions. The basic model, the well-known boson- Josephson 
junction 0, is also analogous to a nonrigid quan- 
tum pendulum. In particular, it was already used for 
observation of the macroscopic quantum tunneling and 
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FIG. 1. (Color online) A schematic representation of the 
setup. Atoms are removed from the left potential well through 
scattering process by an electron or laser beam. 



self-trapping whereas the future proposals include 
also the atomic Mach-Zchnder interferometer [l^l, the 
sensitive weak force detector [l^ and the atomic inter- 
ferometer on a chip [ij]. It was previotisly used also to 
study thermal vs quantum decoherence [15| . 

In our setup, schematically depicted in Fig. [TJ we in- 
troduce an additional element - the controlled removal 
of atoms from one of the two wells of the potential - and 
consider the effect of the dissipation, showing that the 
actual rate of atom losses has a dramatic dependence on 
the dissipation coefficient and on the atomic interactions 
strength. 

We notice that in the earlier experiments [l^ , the Zeno 
effect was observed in a system of cold atoms loaded in 
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an accelerating optical lattice with the magnitude of the 
acceleration varying in time. The macroscopic manifes- 
tation of Zeno and anti-Zeno effects in the Josephson 
junctions subjected to time-dependent perturbations was 
also studied [17[- Also, our statement of the problem, can 
be further developed to describe the experimental setup 
used in [Tsj for the observation of pulsed and continuous 
Zeno effects in an externally driven mixture of two hy- 
perfine states of a ^^Rb BEC, using the destructive mea- 
surement of the population of the states. However, this 
happens in a binary mixture with nonlinear interactions 
between the components, while separation of the two sub- 
systems by a potential barrier in our setting induces only 
a linear coupling between them (i.e. the populations in 
the two wells) due to the quantum tunneling. 

The organization of the article is as follows. In Sec. [IT] 
we deduce the master equation describing our system. 
An exact solution is found in Sec. IIIII for noninteracting 
atoms. The nonlinear Zeno effect is described in Sec. IIVI 
Sec. |V]is devoted to the description of several switching 
regimes which are induced and controlled by the dissi- 
pation. In the Conclusion (Sec. VI) we summarize our 
results and present an outlook. 

II. THE MASTER EQUATION 

We use the simplest boson- Josephson Hamiltonian 
[IS [S HSl describing BEC in an asymmetric double- well 
potential V{r), that is, 

H = - J{a\a2 + alai) + eni + Uin\ + U2n\, (1) 

where aj and (j = 1,2) are boson operators for the 
wave functions Lpi^2 localized at the potential minima, 
rij = ajaj, J is the single atom tunneling rate, e is the 
zero-point energy bias, [/i_2 = 5/2 / d^r<^f 2j where g = 
ATrfi^as/rn, is the s-wave scattering length, and m is 
the atomic mass (for a trap with the two quasidegenerate 
lowest energy levels, Ui ^ U2 = U with a good accuracy). 

The controlled removal of atoms can be realized, e.g., 
by using the experimental technique based on the elec- 
tron microscopy [21I, . A narrow electron beam is di- 
rected to one of the minima of the potential ionizing the 
atoms and the latter are removed from the condensate. 
Alternatively one can use a narrow laser beam. In both 
cases, the interaction with the beam serves as a continu- 
ous measurement tool and can be described in the frame- 
work of the standard Markovian approximation [23| . In- 
troducing the probability p = p{ki, At), where At is the 
time interval and km is a population of the mth well, 
the single-atom removal event can be cast as a quantum 
channel: 

\ki,k2)\0)B.^ Vp\ki - l,k2)\l)B.+ V^~p\ki,k2)\0)R, 

(2) 

where the atoms are removed from well 1, |fci,fc2) = 
^°'/fr' tfr'T' |0) is the ket vector of the BEC state and 



describes the atom counter. Introducing the atom 
removal rate F, we approximate p{ki,At) « TkiAt for 
small At, much less than the tunneling time tqx = ti/ J- 
In particular, in the experiments with the electron mi- 
croscopy, the atom removal rate is computed to be F « 
I(Ttot /e |22| , where / is the current of the electron beam, 
e is the electric charge of the electron, and (Tion is the 
total ionization cross section. 

In terms of the reduced density matrix p, describing the 
condensate alone, the quantum channel ^ is given by 
the Kraus superoperator representation p — > MqpAIq + 
MipMl, where for a small At we have AIq w l — TniAt/2 
and All ~ yTAtai. This leads to the master equation 
in the Lindblad form 

The Lindblad operator !){■) = a(-)a\ — {ni,(-)}/2 has 
the eigenvalues A e {-N,-N + 1/2, -1/2, 0}, where 
N is the total number of atoms. The eigenstates are in 
the product form p = p^^^ (g) p'^^', with p'^^^ corresponding 
to the jth well. Respectively, the right-hand side of Eq. 
([3|), considered as a superoperator, has eigenvalues in the 
form A = lAE — p,, where AE has the range of values of 
the difference between the energy levels of the Hamilto- 
nian while < ^ < FA^. The only stationary state 
(i.e. attractor) is the zero eigenvalue eigenstate, which 
has Pq^-* = |0)(0|, that is, no atoms in the left well. The 
dissipation part is also responsible, besides the removal 
of atoms, for the gradual loss of coherence between the 
wells of the double-well trap, counteracting the effect of 
the quantum tunneling. In its turn, the tunneling de- 
pends on the interactions between the atoms [tI. [sTfloj. 
thus leading to interesting dissipation regimes governed 
by the master equation ([3]). 



III. NONINTERACTING-ATOMS CASE 

Consider first the case oiU = 0, which can be achieved 
by making negligible using the Feshbach resonance. 
Then Eq. ([3]) can be solved explicitly by using the adjoint 
equation (see, Ref. [23j ) 

f4,«,i, + r{44«.-(fi + if)} (4) 

for an observable A [we distinguish an operator solution 
of Eq. dH) with a hat]. First of ah, for aj{t) {j = 1,2) 
Eq. Q is solved by the ansatz aj{t) = Cji{t)ai + Cj2{t)a2 
with the initial condition aj(0) = aj. This gives 

_ A+e^+* - A_e^-* _ A+e^-* - A_e^+* 

J(e^+* - e^-*) 
- 2m 
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where X± 



2h 



The time dependence of the operators is given by 

the Hermitian conjugated expressions. Moreover, the 
time dependence of an arbitrary operator can be found 
in terms of the operators aj (t) and Oj, (t) . This is due to 
the conditional "Leibnitz rule" for the dissipation part of 
Eq. Q: V*[AB] = V*[A]B + AV*[B] valid when either 
[a|, A] = or [ai,B] = 0. Hence, the coefficients Cij{t) 
are sufficient to determine evolution of any observable. 

Setting, for simplicity, e = 0, and assuming that 
initially the condensate is in the ground state — 



(a|+4)"° 



|0) wc obtain 



{N) = 



r 



cos(2m) 



(711 - n^) = -e-2'_ sin(2m)A^o 



(5) 
(6) 



When r < AJ/h, Eqs. ([5]) and (O describe the decaying 
Rabi oscillations with the decay rate r/2. Surprisingly, 
for r ^ AJ/h, the dynamics is characterized by two dif- 
ferent loss rates: the initial stage with the rate r/2 and, 
for times exceeding l/F, a dramatically reduced dissipa- 
tion rate Fqt ~ l^p- This can be explained as follows. 
Consider, for simplicity, the case T 3> J/fi, that is, when 
Tqt < r. We get for t > l/F: 



("i> 



"g-rQTt 



(7) 



Observe that initially the two potential minima are 
equally populated {{rij) = No/2). Equation (O shows 
that there few atoms {Tqt/T ^ 1) in well 1 after the 
time scale t ^ l/F, that is, the system state is close to 
the zero-eigenvalue eigenstate of the dissipation part of 
Eq. ([3]). Further elimination of the atoms in this regime 
occurs via the quantum tunneling from well 2, thus giving 
origin to the dramatically reduced loss rate Fqt- 

The prevention of losses of atoms by a strong dissipa- 
tion resembles the Zeno effect observed experimentally 
in a different setup [31 ■ Specifically, expressing our Tqt 
through the tunneling frequency ujf; = 2J/h we recover 



the decay rate Fqt 



j/F which appears in the con- 



tinuous Zeno effect of Ref. [li 



IV. THE CASE OF INTERACTING ATOMS 

In the nonlinear case one cannot solve Eq. (|1]) exactly. 
The strategy now is to use the mean-field approximation, 
valid in the limit of a large number of atoms. To this end 



we consider the equations for the averaged quantities: 

d{n 



^ = iJ {{a\a2) - (4ai)) - r(ni) 

d{n2) . t// t \ / t \A 
= -iJ [^{a[a2} - (a>i) j 

^ ((ajaa) - (4ai>) = 2i J ((ni) - (713)) 

+2iU ( (nia\a2) + (nia^ai) — {n2a\a2) — (n2ajaij 



(8a) 
(8b) 



-F (^{a\a2) - (4ai>) : 



(8c) 



To obtain a closed system from Eqs. (|8a |) - (|8c|) in the 
nonlinear case one can decouple the fourth-order corre- 
lators as follows: {njO^yajn) « {nj){a^-,ajii) . This proce- 
dure corresponds to the mean-field approximation, that 
is, — > 00, widely used for description of EEC and cold 
atoms for a large number of particles (we have checked 
its validity using the direct quantum Monte-Carlo simu- 
lations). The mean-field variables, z, and q, correspond 
to the averaged quantities: 



{ni) + (712) 

They satisfy the system 
dz 



^J{ni){n2) 



z ^ -, — ^ ; — r, e ^ . . . , g = 



No 



(9) 



^ = -2V1 - z2 sin(/) - ^ (1 - z^) , (10a) 



— = 2^ 



: COS0 + e + 2Aqz, 



(10b) 
(10c) 



Here we use the dimensionless time r = t/tQx and the 
normalized atom removal rate 7 = TIqt = hT/J. The 
parameter A = UN/ J characterizes the atomic interac- 
tions in the condensate 0, Q • 

We have checked, by comparing with the direct 
quantum Monte-Carlo simulations, that solutions of 
Eqs. (|10a|) - pOc|) averaged in the classical phase space give 
an excellent agreement with Eq. ([3]) for iV ^ 100 (whereas 
a good agreement is observed already for N ~ 10). The 
essentially linear dynamics is observed for the interac- 
tion strength A < 1. However, for strong interactions 
(A 1) the self-trapping state [JQ features a dramatic 
reduction of the atom loss rate with increasing interac- 
tion strength (sec Fig. [21 the self-trapping state is in well 
2). The effect docs not depend on the interaction type 
(attractive, A < 0, or repulsive, A > 0) and appears for 
any value of F. 

This effect can be understood as follows. First, one has 
to distinguish between the two types of the self-trapping 
states in the Hamiltonian system: the running-phase 
states with (f>{t) oc t and the fixed-phase states (see Refs. 
d [13, nil). There are the following fixed-phase self- 
trapping states: cos0 — — sgn(A) and z ~ ±Vl — A~2. 
In the dissipative case, z = —1 defines an invariant sub- 



set of system (jlOp representing the mean-field description 
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FIG. 2. The dissipation dynamics of the self-trapping state. 
Here A = 25, e = and 7 = 2. (a) The ratio of the number of 
the remaining atoms. The inset compares the analytical esti- 
mate for the nonlinear dissipation rate (dashed line) with the 
numerical result, (b) The normalized population imbalance 
{z) (the dashed lines show the dispersion), (c) The evolution 
of the phase (</>); the inset gives the phase fluctuations Ae'*". 
The initial state is a Gaussian distribution in the phase space 
with (z) = —0.5 and {(f)) = 0, and dispersions = 0.05 and 
Ae'*^ = 0.1. 



of the zero-eigenvalue eigenstate of the Lindblad term in 
Eq. ([3]), corresponding to the location of all atoms in 
well 2. Accordingly, the fixed-phase self-trapping state 
with 2; < of the Hamiltonian system shows the new 
nonlinear dissipation rate, right after the self-trapping is 
significantly enhanced by the dissipation: z decreases to 
— 1, while the phase becomes equal to an odd (A > 0) 
or even (A < 0) multiple of tt (the phase fluctuations, 
described by the expression Ae*"^ = [1 - |(e''>)|2]V2 
decay to zero). 

Further insight can be gained by visualizing the atomic 
interactions as an effective common potential experienced 
by the condensed atoms. Indeed, in the self-trapping 
state (rii) <C (?t-2); hence, the nonlinear term in the av- 
eraged Hamiltonian ([T]) can be simplified 

U [{nl) + {nl)] « U{n2? « U{N) - 2[/(A^)(ni),(ll) 

where we have used that the fluctuations of z are small 
(see Fig. that is, « {'^2)'^ ■ The total number 

of atoms evolves adiabatically, since the dissipation rate 
is very small. Then, the last term in Eq. dill) , which 
is proportional to (ni), is but a simple renormalization 
of e: £ £ - 2U{N) in Eq. dg [respectively, ((TOb)) ]. 
This allows one to derive the nonlinear dissipation rate 
of the self-trapping state. Considering the symmetric 
trap (e = 0) and setting enl = — 2[/(A^), provided that 
fi^r^ + 4e2jL > J2, we get for t > l/T 



Fnl 



4j2r 



(12) 



This estimate turns out to be in excellent agreement with 
the numerical results (at the initial stage of the decay 
of the self-trapping state), as it is shown in Fig. [2] (if 
one uses in Eq. (jl2p the actual numerical number of 



atoms (N) remaining in the self-trapping state). In the 
regime with Enl 3> hV the nonlinear dissipation rate (jl2p 
reduces to Fnl ~ J^F/e^l. The preceding inequality 
condition also means that Fnl ^ F, that is, the actual 
dissipation rate is dramatically reduced. This inhibition 
of the losses for A S> 1 can be viewed as a nonlinear 
Zeno effect which, in contrast to the usual Zeno effect 
[2^, appears for arbitrary F. 

To conclude this section let us make two remarks on 
the nonlinear decay rate given by Eq. (HU. First, we 
notice that while the analogy between the effect of two- 
body interactions and the bias e in pure linear system 
worked well in the estimate for the decay rate in of the 
nonlinear Zeno effect, we emphasize that the role of these 
factors in the dynamics under consideration is very dif- 
ferent. The linear bias is an external factor which defines 
the constant decay rate of the linear Zeno effect, that is, 
of a condensate of noninteracting atoms. The "induced 
bias" due to the two-body interactions, that is, —2U (N), 
by itself is a dynamical quantity: It slowly changes with 
time, resulting, after all, in the change of the decay rate. 
In other words, after a sufficiently long time, sufficient 
for a significant loss of atoms, the nonlinear decay rate 
will coincide with the linear one. Second, the decay rate 
given in Eq. (jl2p has a broader application than just giv- 
ing the decay of the self-trapping state (thus, it is a new 
nonlinear effect, different from the self-trapping itself). 
For instance, it gives a correction due to the nonlinearity 
to the decay rate of the linear Zeno effect also for a small 
A. Indeed, there are two conditions of validity of Eq. 
P^ : fi^T^ + 4eNL > J^i satisfied also by taking large 
F 3> J/h (and an arbitrary A), and that there few atoms 
in well 1, which is the stage of the Zeno effect. 



QUANTUM SWITCHING INDUCED BY 
DISSIPATION 



Quantum switching can be induced by simply remov- 
ing atoms for a short time in a switch-on manner. For 
instance, application of the atom removal for a time in- 
terval on the order of 1 /F draws the system close to the 
self-trapping state with z « —1. This is illustrated in 
Fig. [31 where about 50% of initially loaded atoms re- 
main in a self-trapping state, whereas until the action of 
dissipation, the system was in the Josephson oscillations 
regime. 

By directing the dissipation tool to the well where the 
self-trapping state is located and by varying the time of 
application of the dissipation and the other parameters, 
for example. A, one can obtain the switching between the 
self-trapping states in the two wells of the double-well 
trap. Fig. |DJa), or induce the switching to the macro- 
scopic quantum tunneling regime. Fig. IH^b). The nona- 
diabatic variation of the dissipation induces large phase 
fluctuations, Ae'"^ = [1 - |(e''^)|2]i/2 ^ 1^ see Fig. [3)^b), 
in contrast to the continuous action of a constant dissi- 
pation, where the phase fluctuations rapidly decay. 
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FIG. 3. Dissipation-induced self-trapping state. Here A = 25 
and an initial Gaussian distribution in the phase space with 
{z) = 0.35, {4>) = 0, Az = 0.03 and Ae"* =0.1 was used. 
In the main panel, the solid line gives the average and the 
dashed lines give the dispersion, 7(f) is given in inset (a) and 
in inset (b) the phase fluctuations are shown. 
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FIG. 4. Dissipation-induced switching. In (a), A = 15 and 
the dissipation acted in the interval At = 0.5tQT leaving 50% 
of atoms in the system. In (b), A = 25 and the dissipation is 
switched for an interval At — O.Sfqx leaving 55% of atoms in 
the system. The dissipation is switched on at t/tgr ~ 1 and 
has 7max = 3. The solid and dashed lines denote the average 
and dispersion, respectively. The initial Gaussian distribution 
has {z) = 0.6, {(j)) = 0, Az = 0.02 and Ae"* = 0.2. 



VI. CONCLUSION 

Selective removal of atoms by a directed external beam 
of particles can be viewed as a continuous measurement 
or, alternatively, as an induced controlled dissipation of 



a mesoscopic system. Combining these, apparently dif- 
ferent interpretations of the interaction between a given 
mesoscopic system and an external system considered as 
a reservoir, with the nonlinearity due to interatomic in- 
teractions, opens remarkable possibilities for observation 
of the fundamental quantum phenomena in open systems 
on one hand, and on the other hand, for preparation and 
controlled manipulation of mesoscopic systems, for ex- 
ample, a Bose-Einstein condensate. More specifically, we 
have shown that while the controlled dissipation atten- 
uates losses of atoms (already known phenomenon, in- 
terpreted also as the celebrated Zeno effect), when com- 
bined with the strong atomic interactions it results in an 
essentially new regime, where the atomic decay rate is 
practically zero, due to emergence of a quasi-self-trapping 
state. Such a state still loses atoms (unlike the case of the 
standard self-trapping state at a fixed number of atoms) 
and in this sense, the nonlinear Zeno effect can be viewed 
as a signature of the self-trapping state. This new effect 
of drastic attenuation of the losses of atoms by the atomic 
interactions, which appears for any value of the dissipa- 
tion parameter, can be viewed as the nonlinear Zeno ef- 
fect. Moreover, application of the dissipation during a 
short interval of time opens new possibilites for control 
over the condensate, for example, inducing the switching 
between two self-trapping states in the two wells of the 
double- well trap or between a self-trapping state and the 
macroscopic quantum tunneling regime. 

We have considered only the simplest nontrivial trap, 
the double-well potential. In general, one should expect 
much richer behavior in the multiwell potentials and, 
moreover, in the case of BEC loaded in the optical lat- 
tices, manipulated by a local dissipation. 

The macroscopic nature of the described phenom- 
ena suggests that they can be observed in any nonlin- 
ear physical system which is described by the nonlinear 
Schrodinger-like equation and has at least two different 
equilibria, one of which is subjected to dissipativc losses, 
that is, in a fairly generic setup. For instance, there is 
a very close analogy with the nonlinear optics of Kerr 
media, in particular, the nonlinear optical fibers (say 
the standard Kerr fibers or hollow-core fibers filled with 
atomic gases), suggesting that the nonlinear Zeno effect 
can be observed also in the realm of nonlinear optics in 
the currently available experimental settings. 
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